Toxicity and Environmental Risk Assessment of Polycarbamate and Its Main Components to Marine Algae and Crustaceans

Polycarbamate is commonly used as an antifoulant coating on fishing nets in Japan. Although its toxicity to freshwater organisms has been reported, its toxicity to marine organisms is currently unknown. We conducted algal growth inhibition and crustacean immobilization tests to assess the effects of polycarbamate on marine organisms. We also evaluated the acute toxicity of the main components of polycarbamate, namely, dimethyldithiocarbamate and ethylenebisdithiocarbamate, to algae, which are the most sensitive tested organisms to polycarbamate. The toxicities of dimethyldithiocarbamate and ethylenebisdithiocarbamate partially explain that of polycarbamate. To assess the primary risk, we derived the predicted no-effect concentration (PNEC) for polycarbamate in a probabilistic manner using species sensitivity distributions. The 72 h no observed effect concentration (NOEC) of polycarbamate to the alga Skeletonema marinoi–dohrnii complex was 0.45 μg/L. The toxicity of dimethyldithiocarbamate may have contributed up to 72% of the toxicity observed for polycarbamate. The fifth percentile of hazardous concentration (HC5) derived from the acute toxicity values was 0.48 μg/L. Comparison of previously reported environmental polycarbamate concentrations in Hiroshima Bay, Japan, to the PNEC estimated using the minimum NOEC and HC5 suggest that polycarbamate currently poses a high ecological risk. Therefore, reducing the risk by restricting polycarbamate use is necessary.


Introduction
Some antifouling compounds have been used on fishing nets and ship hulls as alternatives to organotin antifoulants, which are restricted compounds [1]. Among these antifouling compounds, polycarbamate (PC) is the most commonly used biocide on fishing nets in Japan, accounting for more than 71% of total emissions in the fiscal year 2020 as reported by the Japanese PRTR [2]. This article focuses on the ecological toxicity of the PC that has been widely used on fishing nets over the past two decades.
PC, a dithiocarbamate (DTC) fungicide, was originally developed as a pesticide to control fungal diseases in various field crops [3]. DTCs have long been used for agricultural purposes, and agricultural runoff of these chemicals into surrounding water bodies is inevitable; therefore, their toxicity to freshwater organisms has been evaluated and reported [4][5][6]. However, only a few studies have discussed the toxicity of PC to fish in the marine environment [7,8], and there are no reports on toxicity to algae or crustaceans. Because PC has also been used as an antifouling biocide in coastal waters, its impact on marine organisms is a matter of concern.
Dimethyldithiocarbamate (DMDC) and ethylenebisdithiocarbamate (EBDC) are the main components of PC, considering that PC is formed by two molecules of DMDC and one molecule of EBDC combined with zinc ( Figure 1). PC is hydrolyzed into DMDC and EBDC in aquatic environments because the zinc in PC dissolves in water and the cross-links The species sensitivity distribution (SSD) approach (developed by the United Environmental Protection Agency, Washington, D.C., USA [US EPA]) has been emp to evaluate the ecological risks of pesticide pollution by statistically describing the dependent effects of contaminants at the species level [11]. Using this approach, th ardous concentration that affects 5% of the species (HC5) has been identified as portant indicator of aquatic safety; HC5 (with a 95% confidence interval) is often nized as the threshold for aquatic life [12]. To increase the robustness of the indic has been proposed that the HC5 should be calculated from the model-averaged SS tained from multiple species distributions [13] using an R library [14]. The guidelin SSD approaches recommend using an uncertainty factor to obtain the predicted no concentration (PNEC) from the HC5 values [15][16][17]. The utilization of probabilisti yses to characterize risks has recently become a popular approach and is now the o method for assessing the ecological risks of chemicals in many countries, includi USA, the European Union, Australia, and New Zealand [18][19][20]. However, the risk by PC has not been characterized using probabilistic analysis to date.
The primary aim of this study was to evaluate the acute and chronic toxicities to the marine algae Dunaliella tertiolecta Butcher, Skeletonema marinoi-dohrnii comp complex of Skeletonema marinoi Sarno & Zingone and S. dohrnii Sarno & Kooistra Tetraselmis tetrathele (West) Butcher and the crustaceans Japanese blue crabs (Po trituberculatus Miers) and Tigriopus japonicus Mori. In addition, the acute toxicity of D and EBDC was evaluated using the S. marinoi-dohrnii complex, which is the most se organism to PC out of all tested organisms. The secondary aim of this study was to the ecological risk of PC toxicity. The PNEC was estimated using conventional and abilistic approaches based on the toxicity values derived in the present study. Th dicted environmental concentrations of PC were estimated from the previously re DMDC concentrations in Hiroshima Bay [9], assuming that all detected DMDC mol were derived from PC. The ecological risks from the derived PNEC value and the P The species sensitivity distribution (SSD) approach (developed by the United States Environmental Protection Agency, Washington, D.C., USA [US EPA]) has been employed to evaluate the ecological risks of pesticide pollution by statistically describing the dosedependent effects of contaminants at the species level [11]. Using this approach, the hazardous concentration that affects 5% of the species (HC 5 ) has been identified as an important indicator of aquatic safety; HC 5 (with a 95% confidence interval) is often recognized as the threshold for aquatic life [12]. To increase the robustness of the indicator, it has been proposed that the HC 5 should be calculated from the model-averaged SSD obtained from multiple species distributions [13] using an R library [14]. The guidelines for SSD approaches recommend using an uncertainty factor to obtain the predicted no-effect concentration (PNEC) from the HC 5 values [15][16][17]. The utilization of probabilistic analyses to characterize risks has recently become a popular approach and is now the official method for assessing the ecological risks of chemicals in many countries, including the USA, the European Union, Australia, and New Zealand [18][19][20]. However, the risk posed by PC has not been characterized using probabilistic analysis to date.
The primary aim of this study was to evaluate the acute and chronic toxicities of PC to the marine algae Dunaliella tertiolecta Butcher, Skeletonema marinoi-dohrnii complex (a complex of Skeletonema marinoi Sarno & Zingone and S. dohrnii Sarno & Kooistra), and Tetraselmis tetrathele (West) Butcher and the crustaceans Japanese blue crabs (Portunus trituberculatus Miers) and Tigriopus japonicus Mori. In addition, the acute toxicity of DMDC and EBDC was evaluated using the S. marinoi-dohrnii complex, which is the most sensitive organism to PC out of all tested organisms. The secondary aim of this study was to assess the ecological risk of PC toxicity. The PNEC was estimated using conventional and probabilistic approaches based on the toxicity values derived in the present study. The predicted environmental concentrations of PC were estimated from the previously reported DMDC concentrations in Hiroshima Bay [9], assuming that all detected DMDC molecules were derived from PC. The ecological risks from the derived PNEC value and the PC concentration were investigated.

Polycarbamate (PC)
Actual measured concentrations of PC averaged about 12 ± 3% (mean ± standard error) of the targeted nominal concentrations for the algal growth inhibition tests ( Table 1). The algal growth inhibition tests using D. tertiolecta and the S. marinoi-dohrnii complex satisfied the performance criteria described in the Organization for Economic Cooperation and Development (OECD) test guideline [21], i.e., the growth rate of these algae in the control was higher than 0.92/day during the test period. The 50% effective concentration (EC 50 ) at 72 h and the no observed effect concentration (NOEC) of PC for T. tetrathele were used for reference because the growth rate in the control was lower than 0.92/day during the test period. The 72 h EC 50 of PC for D. tertiolecta, S. marinoi-dohrnii complex, and T. tetrathele were 3.0, 0.72, and 12 µg/L, respectively. The NOEC of PC for D. tertiolecta, S. marinoi-dohrnii complex, and T. tetrathele were 0.48, 0.45, and 1.2 µg/L, respectively (Table 1). The actual measured concentrations of PC averaged about 36 ± 2% of the targeted nominal concentrations for the crustacean immobilization tests (Table 2). In the crustacean immobilization tests of PC, the immobilization in the solvent control group was ≤10% ( Table 2). The 50% lethal concentration (LC 50 ) of PC at 24 h for Japanese blue crab and T. japonicus was 16 and 2.0 µg/L, respectively, under dark conditions at 20 • C.

Main Degradation Products of Polycarbamate (PC)
The actual measured concentration of the sodium salt of DMDC (Na-DMDC), which is a main component of PC, averaged about 15 ± 2% of the targeted nominal concentrations for the algal growth inhibition tests. The 72 h EC 50 and NOEC of Na-DMDC for the S. marinoi-dohrnii complex based on the nominal concentration were 5.2 and 3.5 µg/L, respectively, and the values based on the actual concentration were 0.84 and 0.48 µg/L, respectively (Table 3). In the test of the disodium salt of EBDC (2Na-EBDC), another main component of PC, the chemical was not detected at sufficiently high sensitivity in the test solution. Based on the nominal concentration, the 72 h EC 50 and NOEC of EBDC for the S. marinoi-dohrnii complex were 15 and 10 µg/L, respectively (Table 3).

Species Sensitivity Distribution (SSD)
The SSD approach for PC was carried out using the acute toxicity data (LC 50 and EC 50 ) obtained in the present study and reported data. When two toxicity values were obtained from one species, the geometric mean of the values was used for SSD analysis (Table 4). Acute toxicity values from nine species (three algal species, two crustacean species, and four teleost fish species) were used for deriving SSDs with the cumulative distribution functions using seven distributions, namely, Weibull, log-logistic, log-normal, gamma, log-normal/log-normal mix, log-logistic/log-logistic mix, and Gumbel, and model-fit statistics were calculated (Table 5). A model-averaged SSD was determined based on the Akaike information criterion corrected for sample size (AICc) values using the R package ssdtools [14] and establishing a model-averaged HC 5 (Figure 2). The HC 5 value (0.48 [0.095-5.9] µg/L) was obtained based on the acute toxicity of PC. The results of the Shapiro-Wilk test confirmed the normality of the log-transformed SSD data (p = 0.4976 for the acute toxicity values). The SSD approach was not performed using the chronic toxicity values because at least six toxicity values are required for calculations using the R library, but only five chronic toxicity values (NOEC) were obtained from the present study and published data.    96.0 0.692 0.189 gamma: gamma distribution; lgumbel: Gumbel distribution; llogis: log-logistic distributi llogis_llogis: log-logistic/log-logistic mixture distribution; lnorm: log-normal distributi lnorm_lnorm: log-normal/log-normal mixture distribution; weibull: Weibull distribution.

Toxicity to Marine Algae and Crustaceans
To the best of our knowledge, this is the first study to evaluate the acute and chronic toxicities of PC to marine algae and crustaceans. The acute and chronic toxicities of PC to marine fish have been reported [7,8], and those to marine organisms are shown in Table 4. Algae, especially the diatom S. marinoi-dohrnii complex, are thought to be more sensitive to PC compared to crustaceans and fish, similar to their greater sensitivity to other DTCs [4].
Based on the nominal concentration, the toxicity of PC, Na-DMDC, and 2Na-EBDC to the S. marinoi-dohrnii complex, which is the most sensitive test organism to PC, was 13, 36, and 59 nM, respectively (Tables 1 and 2). PC breaks down into two molecules of DMDC and one molecule of EBDC in water, according to the structural formula of PC ( Figure 1). The contribution of DMDC and EBDC to the toxicity of PC was estimated based on the concept of concentration addition [22], 13 nM (concentration of PC)/13 nM (EC 50 of PC) = 13 nM × 2 molecules (concentration of DMDC)/36 nM (EC 50 of DMDC) + 13 nM (concentration of EBDC)/59 nM (EC 50 of EBDC) = 0.72 + 0.22 = 0.94. Thus, the toxicities of DMDC and EBDC could almost explain that of PC, and the toxicity of DMDC was thought to contribute the most (72%) to the toxicity of PC based on the results of the algal growth inhibition tests using the S. marinoi-dohrnii complex. This study elucidates the toxicity of the degradation products of PC to the marine alga S. marinoi-dohrnii complex. The toxicities of these products to other marine organisms, especially crustaceans and fish, should also be elucidated.

Derivation of Predicted No-Effect Concentration (PNEC)
To the best of our knowledge, this is the first study to estimate the PNECs of PC using probabilistic techniques. Toxicity values from a sufficiently large and ecologically diverse sample of species are needed to estimate the PNEC reliably. In the present study, SSD analyses based on acute toxicity values were carried out using ecotoxicity data obtained from nine species. In addition, the acute toxicity values include the values from species belonging to three trophic levels. The number of toxicity values required to determine SSD varies between guidelines from each nation (e.g., greater than or equal to four through eight) [16][17][18]. Five or more toxicity values are essential for SSD analysis if the ecotoxicity data assumed to follow a log-normal distribution [23] are both cost-effective and highperformance enough to be included in the analysis [24]. At least six toxicity values are required for calculations using the R library ssdtools [14]. Therefore, the HC 5 obtained in the present SSD analysis is considered appropriate for assessing the ecological risk of PC. The derivation of a value corresponding to PNEC from HC 5 varies between the guidelines. HC 5 itself is seldom used as a PNEC. Although an uncertainty factor of 1-5 is applied to HC 5 [15][16][17] in many cases, a default factor of 10 is applied to HC 5 derived from acute toxicity values [16]. Overall, the factor is determined by experts by considering the number of toxicity values and composition of species used in the SSD analysis [15,16]. In the present study, an HC 5 of 0.48 µg/L was obtained from the acute toxicity values. The uncertainty factor of 10 was applied to the HC 5 . A PNEC of 0.048 µg/L was estimated through the probabilistic approach.
The chronic toxicity of PC for marine organisms from two trophic levels was determined, and the lowest toxicity was 0.45 µg/L for the marine alga S. marinoi-dohrnii complex. Conventionally, an uncertainty factor of 50 or 100 is usually applied when chronic toxicity values are obtained from two trophic levels, according to some guidelines, such as the Chemical Substances Control Law in Japan [25] and OECD guidelines [26]. In the present study, 100 was used as the uncertainty factor because this is the most conservative value in the range of recommended values. Thus, the PNEC was calculated to be 0.0045 µg/L. In addition, the toxicity of DMDC to a marine alga was 0.48 µg/L when the same uncertainty factor (100) was adopted, and the reference PNEC was calculated to be 0.0048 µg/L.

Ecological Risk Assessment of Polycarbamate (PC)
In this study, the PNEC of PC for marine organisms was derived using both conventional and probabilistic approaches. Assuming that the detected analyte was derived from PC, our recent investigations estimated the maximum PC concentration in seawater to be 0.11 µg/L in Hiroshima Bay [9], which is higher than the PNECs derived using both the conventional and probabilistic approaches. The maximum DMDC concentration in seawater was reportedly 0.043 µg/L in Hiroshima Bay [9], which is higher than the PNECs obtained using the conventional approach. The ecological risk of PC and DMDC contamination in the Seto Inland Sea is thus currently considered to be high. DTCs including PC are thought to be more stable in seawater than in freshwater because the half-lives of DTCs in solutions increase with increased pH [4]. PCs have been mainly used to coat fishing nets in marine environments [1]. Therefore, the ecological risk of DTCs in the marine environment may be higher than in the freshwater environment.

Scope for Future Research
This study demonstrated that PC is highly toxic to some algae and crustaceans but did not examine the chronic effects of PC on crustaceans or the effects of PC on algae other than microalgae (e.g., macroalgae). The effects of PC on these organisms need to be studied in the future. As DTCs in the environment are broken down into DMDC and EBDC, it is difficult to analyze undegraded DTCs separately, and it is not possible to assess the ecological risk of each DTC without preconditions. Therefore, the toxicity of the degradation products of DTCs such as DMDC and EBDC to marine organisms must be clarified and the ecological risk assessed in more detail.

Organisms
The test organisms used in this study included algae (D. tertiolecta, T. tetrathele, and S. marinoi-dohrnii complex [NIES-324]) and crustaceans (T. japonicus and Japanese blue crab [P. trituberculatus]). Nauplii of T. japonicus (<24 h after hatching) and hatched zoeae of Japanese blue crabs were used in the crustacean immobilization tests.
Dunaliella tertiolecta and T. tetrathele were provided by the Shiogama Branch, Japan Fisheries Research and Education Agency (Shiogama, Japan), and the S. marinoi-dohrnii complex was purchased from the National Institute for Environmental Studies (NIES-324; Tsukuba, Japan). Each algal species was maintained according to the procedure described by Onduka et al. [27] in f/2 medium [28] under static-renewal conditions on a 14:10 h light:dark cycle (light intensity: 39 ± 8 µmol/m 2 /s in the wavelength range 360-670 nm) in a temperature-controlled (20 ± 1 • C) growth chamber (MLR-350; Sanyo, Osaka, Japan). The f/2 growth medium was prepared from natural seawater filtered through sand, activated carbon, and a glass fiber filter (GF/C, Whatman, Maidstone, UK; hereafter referred to as "filtered seawater").
Adult T. japonicus were obtained from the Marine Ecology Research Institute (Tokyo, Japan) and cultured, and T. japonicus nauplii (age < 24 h) from this culture were collected in a 1 L glass beaker containing filtered seawater as described by Onduka et al. [29]. The culture conditions were the same as those used for culturing algae, except that cultured T. tetrathele was added as food. Japanese blue crabs carrying eggs were obtained from the Hamasui Fisheries Company (Hiroshima, Japan). The eggs were placed in glass Petri dishes (90 mm in diameter) containing filtered seawater and incubated at 20 • C overnight with gentle shaking. The hatched zoeae were then used for the toxicity tests.
(special grade, Kyoto, Japan). Other chemicals were obtained from Wako Pure Chemical Industries unless otherwise noted.

Toxicity Tests
Toxicity tests were performed by following the test guidelines of the OECD [21,30] that were modified for marine organisms. Algal growth inhibition tests and crustacean immobilization tests were carried out as described by Onduka et al. [27,31].

Algae
To perform the algal growth inhibition tests, test solutions were prepared by diluting the stock solution (1000 mg/L) and DMSO (solvent control) 3000 times with f/2 medium. The algal growth inhibition tests were carried out in triplicate with a geometric series of nominal PC, Na-DMDC, and 2Na-EBDC concentrations (Tables 1 and 3). The assays were performed using f/2 medium under the same photoperiod and temperature conditions used for the algal cultures and a pH of 8.6 ± 0.6, as described by Onduka et al. [27]. The light source consisted of three ultraviolet-screened fluorescent tubes that yielded a photosynthetic photon flux (360-670 nm) of 41 ± 7 µmol/m 2 /s. During the algal growth inhibition test, the specific growth rate (µ; d −1 ) was calculated as µ = 1/3 × ln (N 72h /N 0h ), where N 0h is the initial cell concentration, and N 72h is the cell concentration after 72 h. The cell concentration was calculated from the relationship between the in vivo fluorescence and the concentration of the alga. The percentage inhibition of growth rate (I) for each treatment was calculated as I = (µ c − µ t )/µ c × 100, where µ c and µ t are the average specific growth rates in the control and exposure groups, respectively. A percentage inhibition of more than 100% means that the cell concentration after 72 h was lower than the initial cell concentration.

Crustaceans
To perform the crustacean immobilization tests, test solutions were prepared by diluting the stock solution (1000 mg/L) and DMSO (solvent control) 2000 times with filtered seawater. The crustacean immobilization tests were performed in a glass weighing bottle (ϕ25 × 25 mm, 6 mL). Five nauplii or zoeae were placed in a well containing 2 mL of filtered seawater. Four wells were prepared for each concentration of the test substances.
The tests were carried out using a geometric series of nominal PC concentrations (Table 2). Animals were not fed during the 24 h test period. The tests were conducted in the dark at 20 • C. Immobilization, which is defined as the inability to swim within 15 s after gentle agitation, was checked at the end of the test. The pH range of the test solutions was 8.0 ± 0.0.

Chemical Analysis
PC and Na-DMDC in the test solutions were analyzed based on the decomposition of PC and Na-DMDC to DMDC followed by methyl derivatization to DMDC-methyl [9]. Water samples were collected and analyzed for the test chemical at 0 and 72 h for the algal growth inhibition test and at 0 and 24 h for the crustacean immobilization tests. The test solutions at the end of the algal growth inhibition test and crustacean immobilization tests were pooled for each exposure group. The test solutions were filtered through a GF/C filter (GE Healthcare, Little Chalfont, UK), and 50 mL (for the algal growth inhibition tests) or 5 mL (for the crustacean immobilization tests) of the solution underwent chemical analysis. The derivatization, extraction, and clean-up methods for PC and Na-DMDC followed the methods described by Hano et al. [9]. Working standard solutions of DMDC-methyl were prepared by diluting the stock solution (100 µg/L in hexane) with hexane. Anthracene-d 10 was used as an internal standard. Elution was performed by adding the internal standard (anthracene-d 10 ; final concentration, 50 µg/L) to each working solution, concentrating them at 40 • C under gentle nitrogen gas to prepare final solutions in 0.5 mL of hexane, and subjecting these solutions to gas chromatography with mass spectrometric detection (GC-MSD; 6890N, 5975 inert; Agilent Technologies, Tokyo, Japan).
For gas chromatographic separation, a 30 m × 0.25 mm (inner diameter) capillary column (HB-5MS; Agilent Technologies) was used. The column was maintained at 60 • C for 2 min, followed by programmed heating to 220 • C at a rate of 20 • C/min. The helium gas flow rate during this procedure was 1.0 mL/min. The mass spectrometer was operated under selected-ion monitoring mode using molecular ions of the target test chemicals (electron impact at 70 eV, ion source temperature, 230 • C). The interface temperature was 275 • C. The ions (m/z) monitored in this mode were DMDC-methyl (88, 135) and anthracene-d 10 (188).
PC concentrations were recalculated from the DMDC-methyl concentrations. For PC analysis, the concentration of PC (p) was calculated as p = d × 0.89/0.41, where d is the concentration of DMDC-methyl [9]. The detection limit (DL) was determined as DL = 10 × signal-to-noise ratio. The DL of DMDC-methyl and PC for 50 mL of the algal test solution were 0.02 and 0.04 µg/L, respectively. The detection limit of the PC for 5 mL of the crustacean test solution was 0.2 µg/L. The geometric mean of the toxicant concentrations was used to estimate the EC 50 , lowest observed effect concentration (LOEC), and NOEC.

Data Analysis
The 24 h and 72 h EC 50 were estimated, and multiple comparison tests were carried out using R software [32] (version 4.2.1), primarily with the drc and multcomp (Dunnett's test) packages [33,34]. Differences between groups were considered significant at p < 0.05. The normality of the log-transformed SSD data was confirmed using R software with the dplyr (Shapiro-Wilk test) package [35].
To develop the SSD curves and estimate the HC 5 values for PC, the maximum likelihood estimation (MLE): model-averaging approach was used [13,36]. This technique is robust because more than one curve can fit the data equally well in some instances, and the choice of the curve impacts the calculation of the HC 5 . By using the model-averaging approach, bias caused by the selection of one curve over another can be avoided. Moreover, the resulting weighted distribution is relatively stable and not greatly influenced by small changes in the dataset [13,36].
The methods and recommendations outlined by Schwarz and Tillmanns [13] and R [32] (version 4.2.1) with R Studio [37] (version 2022.12.0) and the R package ssdtools [14] (version 1.0.2) were used to implement the following cumulative distribution functions based on their ability and robustness of the fitted SSD data: Weibull, log-logistic (llogis), log-normal (lnorm), gamma, log-logistic/log-logistic mix (lnorm_lnorm) log-normal/lognormal mix (lnorm_lnorm), and Gumbel (lgumbel). The quality of the fitted distributions using MLE as the regression method was then assessed for each scenario using the AICc, which is a measure of the relative quality of fit to the dataset, for each distribution [13,38,39]. A model-averaged SSD was determined based on AICc values of multiple distributions using the R package ssdtools [14], and a model-averaged HC 5 was estimated using the model-averaged SSD [13,36]. Parametric bootstrap resampling (using 10,000 iterations) was then performed to calculate confidence intervals across the entire SSD curve, including the HC 5 [13,36].

Conclusions
In this study, the acute and chronic toxicities of PC to marine organisms such as algae and crustaceans were evaluated under laboratory conditions. It was determined that algae, especially the diatom S. marinoi-dohrnii complex, are more sensitive to PC than crustaceans and fish are. In addition, we evaluated the acute toxicity of DMDC and EBDC, which are the main degradation products of PC, to the S. marinoi-dohrnii complex, the most sensitive to PC out of all tested organisms, and demonstrated that the toxicity of DMDC is similar to that of PC and plays a large role in the toxicity of PC. The PNEC of PC for marine organisms was derived through probabilistic and conventional approaches. With regards to the ecological risk of PC, the highest concentration of PC in Hiroshima Bay was 0.11 µg/L, which is higher than the PNECs derived through both the conventional and probabilistic approaches. Therefore, the ecological risk of PC contamination in Hiroshima Bay is currently considered to be high. Hence, stakeholders should formulate strategies for minimizing the risks posed by DTCs including PC, such as, for example, restricted use of DTCs in the bay and the basin of the river running into the bay.

Data Availability Statement:
The data generated in this study are available within the article.